First evidence for molecular interfaces between outflows and 
ambient cloud in high-mass star forming regions? 
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We present new observations of the CepA-East region of massive star formation and describe 
an extended and dynamically distinct feature not previously recognised. This feature is present in 
emission from H2CS, OCS, CH3OH, and HDO at -5.5 km s~^, but is not traced by conventional 
tracers of star forming regions H2S, SO2, SO, CS. The feature is extended up to at least 0.1 pc. 
We show that the feature is neither a hot core nor a shocked outflow. However, the chemistry of 
the feature is consistent with predictions of a model of an eroding interface between a fast wind 
and a dense core; mixing between the two media occurs in the interface on a timescale of 10-50 
years. If these observations are confirmed by detailed maps and by detections in species also 
predicted to be abundant (e.g. HCO"*", H2CO, and NH3) this feature would be the first detection 
of such an interface in regions of massive star formation. An important implication of the model 
is that a significant reservoir of sulfur in grain mantles is required to be in the form of OCS. 

Subject headings: ISM: individual (CepA-East) — ISM: molecules — radio lines: ISM — star formation 
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ABSTRACT 



To explore the interactions of young stellar 
objects (YSOs) with their environments, we re- 
cently carried out mm-wavelength molecular line 
observations towards star forming regions (SFRs) 
with well defined and bright high-velocity compo- 
nents. In particular, we mapped the well known 
Cepheus A (CepA; Garay et al. 1996, and ref- 
erences therein) SFR in several shock-chemistry 
tracers such as H2S, SO2, and HDO (Codella et 
al. 2003, 2005). Those results show that the 



ducing a hot core (Martin-Pintado et al. 2005), 
are also associated with multiple mass loss pro- 
cesses. In particular, beside the already known 
three flows pointing in the SW, NE, and SE direc- 
tions, a fourth outflow flowing towards the South 
has been detected thanks to the shock-chemistry 
tracers. CepA-East can thus be considered an 
ideal laboratory in which to study how outflow 
motions affect the gas, from both the kinematical 
and chemical points of view. Codella et al. (2005) 
have already presented a multi-species and multi- 
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line nim-survcy of the central region of CepA-East 
where the YSOs are located. Using the 30-m 
IRAM antenna, the authors detected emission in 
different transitions of 21 molecular species trac- 
ing a wide range of physical conditions. Analysis 
of these spectra shows that different molecules ex- 
hibit different spectral behaviours and that three 
classes can be distinguished: (i) hot core molecules 
(e.g. HC^*0+, CH3C2H) emitting only at the ve- 
locity of the hot core (-10.7 km s~^) and with no 
line wings, (ii) outflow molecules (e.g. CS, SiO, 
H2S, SO2, and SO) spanning the whole range of 
observed outflowing velocities so that bright wings 
are added to the hot core emission, and (iii) four 
species (OCS, H2CS, HDO, and CH3OH) which 
are associated with wings and which, in addition, 
clearly show a redshifted spectral peak at -5.5 
km s~^, well separated from the hot core peak. 
While the peak at 10.7 km s~^ is tracing the 
high-density material hosting the YSOs and the 
wings are tracing the multiple outflows, the origin 
of the redshifted spectral peak is unknown. The 
Codella et al. (2005) data did not allow us to 
clarify the spatial distribution of this spectral peak 
and to establish if it is tracing a small structure 
or it is related with an extended component. It is 
worth noting that, as far as we know, this is the 
first study to reveal duality in the line-wing pro- 
files observed in outflows driven by YSOs, i.e that 
OCS, H2CS, HDO, and CH3OH (hereafter called 
double-peaked species) have a different behaviour 
with respect to CS, SiO, H2S, SO2, and SO (here- 
after called single-peaked species). This suggests 
that the redshifted spectral peak could be trac- 
ing a different gas component with respect to the 
gas outflowing at the other velocities as well as to 
the hot core gas, and indicates that high angular 
resolution observations are required for a detailed 
analysis. 

In this Letter we present observations which 
clarify the spatial distribution of the redshifted 
spectral peak at -5.5 km s~^. In addition, we 
produce chemical models suggesting that we are 
probably observing the first direct evidence of tur- 
bulent interfaces (i) where the outflow gas im- 
pinges on and detach dense gas, and (ii) in which 
dynamical mixing and diffusion are occurring be- 
tween the mass ejected from a newly- formed mas- 
sive YSO and the interstellar cloud from which 
it was formed. In the following, the kinematical 



component at -5.5 km s ^ will be referred to as 
the I-feature. 

2. The spatial distribution of the I-feature 

The main goal of the observations was to de- 
termine the spatial distribution of the I-feature 
observed towards CepA-East. In order to select 
the best molecular tracers, following Codella et al. 
(2005), we noted that a good compromise between 
high angular resolution, a simple spectral pattern, 
and an intense line emission was represented by 
H2CS, and in particular by its 6i6-5i5 transition 
at 202923.55 MHz {E^^A7 K). Thus, we mapped a 
region of - l'x2' in H2CS(6i6-5i5) on 2004 June 
with the 30-m IRAM radiotelescope of Pico Veleta 
(Spain). We used a sampling of 10" around the co- 
ordinates of HW2 YSO, which is thought to be 
among the main drivers of the CepA-East out- 
flows and it represents the center of the present 
map, whereas a spacing of 20" was chosen to scan 
coarsely the outer part of the cloud. The sys- 
tem temperature, Tgys, was ^ 500 K, the HPBW 
was 12", while the pointing was accurate to within 
2"-3". As spectrometer, we used an autocorrela- 
tor with a configuration providing a velocity res- 
olution of 0.06 km s^^, successively smoothed to 
0.92 km s~^. The spectra were calibrated with 
the standard chopper wheel method and reported 
here in units of main-beam brightness temperature 
{Tmb)'- the average r.m.s. is ^ 20 mK. 

Figure 1 a reports the velocity channel maps of 
the H2CS(6i6-5i5) emission. Each panel shows 
the emission integrated over a velocity interval of 
2 km s~^; the ambient velocity (Vlsr) is -10.7 
km s~^ (e.g. Codella et al. 2005). In summary, 
H2CS(6i6-5i5) emission is associated (i) with the 
central position, where the YSO cluster is located 
and where a hot core has been detected (Martin- 
Pintado et al. 2005); and (ii) with the four outflow 
directions, NE, SE, S, and SW, confirming H2CS 
as a tracer of high-temperature and/or shocked 
regions. In particular, the new H2CS maps con- 
firm that the spatial distribution of the I-feature 
is not limited to the central position tracing also 
the SW and southern outflows, as clearly shown 
by the -5.5 and -3.5 km s~^ panels in Fig. la. 

Examples of H2CS(6i6-5i5) line proflles are 
shown in Fig. lb which compares the spectra ob- 
served at the (0",+10") and (0",-10") offsets with 
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that observed at the (0",0") position (CodcUa et 
al. 2005). Given the HPBW, the three spec- 
tra are sampHng different regions of CepA-East. 
The I-fcaturc is still present in the southern po- 
sition and it is even more redshifted, approach- 
ing the -5 km velocity, thus suggesting the 
presence of a velocity gradient. In order to make 
the same comparison with another tracer of the 
I-feature, we observed the CH30H(5k-4k) tran- 
sition at 241.8 GHz towards the (0",+lG") and 
(0",-10") positions. In this case, Tgys — 650 K, 
the HPBW = 10", and the resulting spectra were 
smoothed to 0.39 km s~^, with a r.m.s. ^ 80 mK. 
Fig. 1 c compares the profiles due to the high exci- 
tation (£;u=115 K) CH30H(54-44 A±) transition 
at 241806.51 MHz, which according to the (0",0") 
spectrum, can be considered as the best tracer of 
the I-feature among the lines of the CH30H(5k- 
4k) pattern. Also in this case, the I-fcaturc is still 
present in the southern position and again it seems 
more redshifted, in agreement with the H2CS data. 
It is worth noting that the two continuous lines 
across the methanol (0",0") spectrum refer to the 
two velocity components of the CH30H(5_4-4_4 
E) transition {E^=122 K): the 5.5 km s^^ (5_4- 
4_4 E) emission is blended with the -10.7 km s~^ 
(54-44 A='=) feature, producing a broad line. Fi- 
nally, Fig. Id reports the HDO(lio In) (£u=47 
K) profiles that were also reported in the Codella 
et al. (2005) paper. Since the angular resolu- 
tion is definitely higher (HPBW=31") wc show the 
spectra taken in the {0" ,+20") and (0",-20") posi- 
tions to sample different parts of CepA-East. The 
HDQ data are well in agreement with the picture 
given above, showing how bright is the I-feature 
as observed towards the southern position. 

We conclude that there is a previously unde- 
tected structure extending toward the South from 
the central position giving rise to emission at - 
5.5 km s-i in lines of H2CS, OCS, HDO, and 
CH3OH. The structure appears to be at least 0.1 
pc in length. 

3. The origin of the I-feature 

In this Section wc consider possible origins for 
the extended I-feature at -5.5 km s^^. Since 
CepA-East is a YSO cluster the possibility that 
this emission comes from a second hot core seems, 
at first sight, likely. However, this cannot be the 



case because: (i) the emission is extended and 
not compact as it would be if it came from a 
hot core (see Fig. 1). In addition, (ii) Figure 2 
shows that what makes the I-fcat\irc pcciiliar is 
that some line ratios such as OCS/SO, H2CS/SO2 
or CH3OH/H2S are definitely much higher than 
those measured towards the 10.7 km compo- 
nent. In other words, although we cannot exlude 
that conventional hot core tracers, including, SO, 
SO2, and H2S, can have some weak emission at - 
5.5 km s~^, they are not tracing the I-feature. Fi- 
nally, (iii) we find that the column densities of the 
double peaked species in the I-feature (CH3OH: 
3 10^5 cm-2; OCS: 5 lO^^; H2CS: 1 lO"; HDO: 
6 10^*; Codella et al. 2005), calculated assuming 
extended emission, arc lower that those typically 
found in hot cores (Jacq et al. 1990; Millar & 
Hatchell 1998). We conclude that there are both 
morphological and chemical reasons for excluding 
a hot core as the origin of the I-feature. 

We have also investigated the possibility that 
the I-feature is simply the emission from a shocked 
molecular outflow. This possibility has been pro- 
posed in the Codella et al. (2005) paper, where 
the observations were compared with the theoret- 
ical calculations recently reported by Wakelam et 
al. (2004). However, a closer inspection of such 
model suggests that the theoretical abundanc;cs of 
the single-peaked species are not reproduced and 
that both single and double peaked species abun- 
dances, once converted into column densities (i.e. 
taking into consideration the geometry) are more 
typical of hot cores densities and hence too high. 
In addition, it is difficult to understand why stan- 
dard shock tracers like SiO. SO2, and H2S do not 
show the same profile observed with the double- 
peaked species, as suggested by shock-chemistry 
models and as happens in the well-studied chemi- 
cal rich outflows BHR71 and L1157 (Garay et al. 
1998; Bachiller & Perez Gutierrez 1999). 

A new explanation is therefore required. We 
now investigate the possibility that the I-feature 
arises from a molecular interface between the out- 
flow gas and the ambient medium. Highly tur- 
bulent interfaces should be present in all environ- 
ments where jets or outflows interact with the sur- 
rounding molecular clouds. Such interfaces have 
been shown to be characterized by a chemistry 
which is quite distinct in nature from that of typ- 
ical dense cores and would not be achieved by 



3 



any modification of conventional cold cloud chem- 
istry. All chemical interface models, whether with 
turbulent mixing or diffusion, explore the conse- 
quences of mixing warm largely ionized gas with 
cold dense and mainly neutral gas (e.g. Lim, Rawl- 
ings, & Williams 2001; Rawlings & Hartquist 
1997). There is as yet no direct evidence of the 
existence of such interfaces in high-mass SFRs, 
while in low-mass SFRs, observations of enhanced 
HCO+ towards several Class objects can be in- 
terpreted as coming from the walls of the outflow 
cavity (Hogerheijde et al. 1998; Rawlings, Tay- 
lor, & Wilhams 2000; Redman ct al. 2004). 
In order to test this hypothesis we make use of 
the Viti et al. (2002) interface model and at- 
tempt to reproduce the different trends displayed 
by the single-peaked and double-peaked species. 
Note that since the model does not include the 
deuterium chemistry, we could not calculate the 
abundance of HDO. 

The main physical characteristics often adopted 
for an interface are: low visual extinctions (< 1.5 
mags); high gas densities (riHa — 10**-10^ cm~^), 
high radiation fields compared to the mean in- 
terstellar radiation field (here called one Habing), 
short lifetimes (< 100 yr) and relatively high tem- 
peratures (Tkin > 100 K) due to non-dissociative 
shocks associated with the highly turbulent inter- 
face between the outflow and the molecular core. 
Note that although timescales may be very short, 
the continual erosion of dense material by the wind 
or jet re-supplies the interface, so a near steady- 
state pertains. The temperature and density of 
the region, as deduced from the molecular obser- 
vations (Codella et al. 2005), are ~ 200 K and 
lO^-lO'^ cm~^, respectively. We have therefore 
run a small grid of models using these values and 
varying the visual extinction (0.5-1.5 mags) and 
the radiation fields (5-100 Habing). A further pa- 
rameter that we investigate is the form that sulfur 
takes once depleted on the grains, before the man- 
tles are evaporated. We have considered 100% in 
solid H2S, 100% in solid OCS and several mixtures 
in between. Note that while it is often assumed 
that most simple species, once frozen out, react 
with hydrogen atoms to form saturated species, 
recently doubts have been raised regarding the 
form that sulfur takes once it depletes onto the 
grains. In particular, the detection of solid OCS 
(Palumbo, Geballe, & Tielens 1997) and the non- 



detection of H2S ices (van Dishoeck & Blake 1998) 
have raised the possibility that OCS may be the 
main reservoir of sulfur on the grains. It is not un- 
realistic to assume that before evaporation most of 
the sulfur is in the form of OCS: studies of the ef- 
fects of UV on grains (e.g. Mannella, Palumbo, & 
Baratta 1997) have shown that species such as 
H2O for example, are easily dissociated by radi- 
ation; once the oxygen is freed, it is believed to 
react with CO, which is the third most abundant 
species on the grains, to form CO2. It is not un- 
reasonable to assume that a similar process occurs 
with H2S: i.e. that once atomic sulfur is freed it 
reacts efficiently with CO to form OCS. 

The aim of this exercise is not to match per- 
fectly the observed column densities but to see 
whether there is an epoch when a molecular 
outflow-ambient interface could give rise to high 
abundances of CH3OH, OCS and H2CS while SO, 
SO2 and H2S remained low. Figure 3 shows one 
of our best fit models, where, in general, the be- 
haviour is as expected: the double-peaked species 
show a clear anti-correlation with respect to the 
single-peaked family. Note that in any case the 
moleciilar abundances of all the species ultimately 
(> 100 yr) decline as the radiation field induces 
photodissociation, but, as mentioned earlier, the 
interface is turbulent and therefore a continu- 
ous erosion of material will constantly replenish 
the interface gas. A radiation field larger than 
10 Habing seems to destroy most of the doubly- 
peaked species too early. A density lower than 10^ 
cm~^ does not produce enough methanol. But the 
most interesting constraint is the amount of sulfur 
(at least 90%) that needs to be in solid OCS before 
evaporation occurs. Thus, these calculations sug- 
gest that the observations are consistent with (but 
do not prove) a model of an outflow- ambient tur- 
bulent interface with the following physical condi- 
tions: - 10^ 10^ cm-'^ Tkin - 200 K, and a 
radiation field between 5 and 10 Habing. Accord- 
ing to this modelling, our observational data may 
be evidence of the presence; of a molecular inter- 
face in high mass star forming regions. As well as 
indicating high abundance of the double-peaked 
species, our model: 

1. requires high abundances of solid OCS, and 

an entrainment of sufficient material into the 
interface on timescales of 10-50 years. While 
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it is feasible' that OCS on ices is abundant 
in an environment where the UV radiation 
field is strong, it is clear that experimental 
as well as observational studies on ices are 
desirable to confirm such picture; 

2. predicts high abundances of other species, 
including HCO+, H2CO, and NH3. To con- 
firm the presence of a molecular interface be- 
tween outflow and ambient gas, we thus sug- 
gest further observations to map CepA-East 
in these species as well as OCS lines. 
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Fig. 1. — (a) Channel map of the H2CS(6i6-5i5) emission towards CepA East. Each panel shows the emission 
integrated over a velocity interval of 2 km s~^ centred at the value given in the bottom left corner. The 
thick boxes point out the main velocities of the I- feature (sec text). The position offsets are relative to Q!2ooo 
= 22^ 56™ 17^9, (^2000 = +62° 01' 49'.'7, which are the HW2 coordinates. The empty circle shows the IRAM 
30-m beam (HPBW; 12"), while the small crosses mark the observed positions. The triangles stand for the 
VLA 2 cm continuum components which trace two strings of sources arising in shocks (Garay ct al. 1996). 
The black arrow points out the direction of the 0.6 pc coUimated H2S outflow discovered by Codella et al. 
(2003). The contours range from 80 to 1120 mK km s~^. The first contour and the step correspond to about 
4 a (where a is the r.m.s. of the map), (b) and (c): H2CS(6i6 -815) {E^=47 K) and CH30H(54-44 A±) 
{£^=115 K) spectra in Tmb scale observed at the (0",+10") and (0",-10") offsets and compared with those 
observed at the (0",0") position (Codella et al. 2005). The HPBWs of the CH3OH and H2CS observations 
are 10" and 12", respectively. The H2CS spectra have been truncated to stress the I-featurc, pointed out by 
a dashed line as well as the ambient velocity, -10.7 km s~^. The two continuous line in the (0",0") CH3OH 
panel refer to the two velocity components of the (5_4-4_4 E) transition (£Ju=122 K). (d): Same as (b) and 
(c) for the HDO(lio In) (£u=47 K) emission (Codella et al. 2005). Since in this case the HPBW is 31", 
the spectra taken at the (0",+20") and (0",-20") are shown. 
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Fig. 2. — Molecular line profiles observed towards CepA East (Codella et al. 2005): species and transitions 
are reported. Two different kinds of line profiles are observed: (i) H2S, SO2, SO, and CS (lower panels) 
are associated with extended wings spanning the whole range of outflow velocities, whereas (ii) H2CS, OCS, 
HDO, and CH3OH (upper panels) are associated with wings and in addition show a redshifted secondary 
peak at -5.5 km s^^ (called I-fcaturc in the text), well separated from the ambient velocity ( 10.7 km s""'^). 
The dashed lines stand for these two velocities. The spectra in the lower panels have been truncated in order 
to focus the attention on the spectral wings and to stress the absence of the secondary peak clearly observed 
thanks to the H2CS, OCS, HDO, and CH3OH emissions. 
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rijj =10^ cm ^ , 5 Habing 
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Fig. 3. — Column densities of the double-peaked (H2CS, OCS, CH3OH) and the single-peaked species (H2S, 
SO2, SO, CS) as a function of time for one of the best fit models of molecular interfaces (see Sect. 3): the 
gas density is 10^ cm~^, while the radiation field is 5 Habing. Note that the column density of H2S is less 
than 10^^ cm~^ and hence it is undetectable. Before evaporation, 100% of the sulfur on the grains is in the 
form of OCS. 
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